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(57) Abstract: ABSTRACT The invention relates to apparatus, methods and computer program code for equalisation. A 
soft-in-soft-out (SISO) equaliser for use in a receiver for receiving data from a transmitter configure to transmit data from a plurality 
of transmit antennas simultaneously is described. The equaliser comprises at least one received signal input for inputting a received 
signal; a plurality of likelihood value inputs, one for each transmit antenna, for inputting a plurality of decoded signal likelihood 
values from a SISO decoder, a processor configure to determine from the plurality of signal likelihood values an estimated mean 
and covariance value for a signal from each of the transmit antennas; an expected signal determiner coupled to the processor to 
detennine an expected received signal value using the mean values; a subtracter coupled to the received signal input to subtract the 
expected received signal value from the received signal to provide a compensated signal; a filter coupled to the sub tractor to filter the 
compensated signal to provide a plurality of estimated transmitted signal values, one for each the transmit antenna; a filter coefficient 
determiner coupled to the processor to detennine coefficients of the filter using the covariance values; and an output stage coupled to 
the filter to output a plurality of transmitted signal likelihood values, one for each the transmit antenna, for input to the SISO decoder. 
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1 

DESCRIPTION 
EQUALISATION APPARATUS AND METHODS 

5 Technical Field 

This invention relates to apparatus, methods and computer program code 
for equalisation in communication systems in which a receiver receives signals 
from a plurality of transmit antennas, and in particular where the plurality of 
transmit antennas is associated with a single transmitter. 

10 Background Art 

A particular problem arises in a communications links where a transmitter 
with more than one transmit antenna is employed since signals received from 
different transmit antennas interfere with one another. This results in so-called 
multi-stream interference (MSI) and causes decoding difficulties. The potential 

1 5 advantage, however, is greatly increased throughput (that is, a higher bit rate) for 

such a communications link. In this type of MIMO (Multiple-input Multiple- 
output) communication link the "input" (to a matrix channel) is provided by the 
transmitter's plurality of transmit antennas and the "output" (from a matrix 
channel) is provided by a plurality of receive antennas. Thus each receive antenna 

20 receives a combination of signals from all the transmitter's transmit antennas 
which must be unscrambled. 

A typical wireless network comprises a plurality of mobile terminals (MT) 
each in radio communication with an access point (AP) or base station of the 
network. The access points are also in communication with a central controller 

25 (CC) which in turn may have a link to other networks, for example a fixed 

Ethernet-type network. Until recently considerable effort was put into designing 
systems so as to mitigate for the perceived detrimental effects of multipath 
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propagation, especially prevalent in wireless LAN (local area network) and other 
mobile communications environments. However the described work G.J. Foschini 
and M J. Gans, "On limits of wireless communications in a fading environment 
when using multiple antennas" Wireless Personal Communications vol. 6, no. 3, 

5 pp. 31 1-335, 1998 has shown that by utilising multiple antenna architectures at 
both the transmitter and receiver (a so-called multiple-input multiple-output 
(MIMO) architecture) vastly increased channel capacities are possible. Attention 
has also turned to the adoption of space-time coding techniques (in OFDM, space- 
frequency coding) for wideband channels. Typically channel state information 

10 (CSI) for maximum likelihood detection of such coding is acquired via training 
sequences and the resulting CSI estimates are then fed to a Viterbi decoder. 

A technique for space-time code detection in a MIMO system based upon 
the use of periodic pilot sequences and interpolation filters is described in A. 
Naguib, V. Tarokh, N Seshadri and A. Calderbank "A space-time coding based 

15 model for high data rate wireless communications" IEEE J-SAC vol. 16, pp. 1459- 

1478. Oct 1998. 

FIG. 1 shows a MIMO communication system 100 in the context of which 
the technique described in Naguib et al operates. An information source 101 
provides an information symbol s(l) at time 1 to a space-time encoder 102 which 

20 encodes the symbol as N code symbols Ci(l) C2O). . . cn(1\ each of which is 

transmitted simultaneously from one of transmit antennas 104. A plurality M of 
receive antennas 106 receives respectively signals ri(l), . . . sm{\) which are input to 
receiver 108. The receiver 108 provides on output 1 10 an estimate s(l) of the 
encoded transmitted symbol s (I). There is a plurality of channels between the 

25 transmit and receive antennas, for example all channels with two transmit antennas 

and two receive antennas. The technique described in Naguib et al uses periodic 
pilot sequences in the transmitted signal to estimate the time varying responses of 
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these channels. 

Third generation mobile phone networks use CDMA (Code Division 
Multiple Access) spread spectrum signals for communicating across the radio 
interface between a mobile station and a base station. These 3G networks are 

5 encompassed by the International Mobile Telecommunications IMT-2000 standard 
( www.ituintV Collectively the radio access portion of a 3G network is referred to 
as UTRAN (Universal Terrestrial Radio Access Network) and a network 
comprising UTRAN access networks is known as a UMTS (Universal Mobile 
Telecommunications System) network. The UMTS system is the subject of 

10 standards produced by the Third Generation Partnership Project (3GPP, 3GPP2), 
technical specifications for which can be found at www.3gpp.org . Fourth 
generation mobile phone networks, although not yet defined, may employ MIMO- 
based techniques. 

In practical data communication systems multipath within a channel results 
15 in intersymbol interference (ISI), which is often corrected with a combination of 
equalisation and forward error coding. For example a linear equaliser effectively 
convolves the received data with an inverse of the channel impulse response to 
produce data estimates with ISI substantially removed. An optimal equaliser may 
employ maximum likelihood (ML) sequence estimation or maximum a priori 
20 estimation (MAP), for example using a Viterbi algorithm. Where data has been 

protected with a convolutional code a soft input Viterbi decoder may be employed, 
usually together with data interleaving to reduce the effects of burst errors. Such 
approaches provide optimal equalisation but become impractical as the symbol 
alphabet size and sequence length (or equivalently channel impulse response 
25 length) increases. 

Turbo equalisation achieves results which are close to optimal but with 
substantially reduced complexity compared to non-iterative joint channel 
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equalisation and decoding. Broadly speaking turbo equalisation refers to an 
iterative process in which soft (likelihood) information is exchanged between an 
equaliser and a decoder until a consensus is reached. The effect of the channel 
response on the data symbols is treated similarly to an error correction code and 

5 typically a soft output Viterbi algorithm (SOVA) is used for both. Again, however, 

such techniques are impractically complex for large delay spreads and symbol 
alphabets, particularly as several processing iterations may be needed to achieve 
convergence for a single data block. These difficulties are significantly 
exacerbated where signals from more than one transmit antenna must be 

10 disentangled and equalised, with a different channel response for each transmit 

antenna or transmit-receive antenna pair. 

A paper by Tuchler et al. (Minimum Mean Squared Error Equalization 
Using A-priori Information, Michael Tuchler, Andrew Singer, Ralf Koetter, IEEE 
Transactions on Signal Processing, vol. 50, pp. 673 - 683, March 2002) describes 

1 5 a simplified approach to turbo-equalisation where a single transmit antenna is 
employed. In this paper the conventional MAP equaliser is replaced by a linear 
equaliser (that is by a linear or transversal filter) with filter coefficients which are 
updated using a minimum mean square error (MMSE) criterion evaluated over 
both the distribution of noise and the distribution of symbols. A linear estimate x n 

20 of a transmitted symbol x n is determined using an observation z n via the equation 
i„= a„ H z n + b n where superscript H denotes the Hermitian operator and a n and b n 
are the coefficients of the estimator (strictly the estimate should be termed affine 
rather than linear because of the constant b n ). The coefficients are chosen to 
minimise the MSE cost E(\x n - x„ | 2 ) where £(•) denotes a mean or expectation 

25 value. As information is fed back to the equaliser from the error correction 

decoder the filter coefficients change with time and are thus recomputed for each 
data symbol to be estimated. A related technique is described in WO 02/15459. 
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5 

The contents of the Tuchler et al. paper are helpful for understanding the 
present invention, which builds upon and extends this work and, in particular, 
sections II and III of this paper are specifically hereby incorporated by reference. 

Another paper, Tetsushi Abe and Tad Matsumoto, "Space-Time Turbo 
5 Equalization and Symbol Detection in Frequency Selective MEMO Channels" in: 
Proc. Veh. Techn. Conference, ffiEE VTS 5 th Vol 2. pg 1230-1234, 2001 
describes the application of turbo equalisation to a MIMO system with a plurality 
of users (transmitters). However the simplified approach described in this paper is 
only suitable for BPSK and not for other (higher) modulation schemes used, for 
10 example, by wireless LAN networks. Therefore improved techniques, and in 

particular improved algorithms for systems with multiple antenna transmitters are 
still desirable. 

Other background prior art may be found in US 2002/01 10188 which also 
employs adaptive linear filtering in SISO (soft-in soft-out) turbo equalisation, in 

15 Iterative (Turbo) Soft Interference Cancellation and Decoding for Coded CDMA, 

Xiaodong Wang, H. Vincent Poor, IEEE Trans. Comms. Vol 47, No. 7, July 1999, 
pp. 1046-1061, and in Iterative Receivers for Multiuser Space-Time Coding 
Systems, Ben Lu and Xiaodong Wang, IEEE Journal On Selected Areas in 
Communications, Vol. 18 No. 11, Nov 2000. pp. 2322-2335. In these documents 

20 the references to "soft" information are to information relating to the likelihood of 
a particular bit, symbol or signal value as opposed to, for example, hard 
information resulting from a hard bit decision defining a bit as a logic one or zero. 
More complex trellis-based and Bayesian-based SISO equaliser components have 
also been utilized for iterative turbo-equalization where the equalizer and decoder 

25 iteratively exchange transmitted symbol likelihood information to improve the bit- 

error-rate of the receiver. 

Filter-based equalisers using MMSE or zero-forcing criteria, 
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Bayesian-based equalisers and trellis-based equalisers such as (MAP) and 
Maximum Likelihood Sequence Estimation (MLSE) equalisers have all been 
employed for wideband MIMO systems to mitigate both ISI and MSI. It is also 
known to use iterative turbo equalization technology for equalisation in MIMO 

5 systems. However turbo equalisers employing trellis or Bayesian-based SISO 

equalisation have a high computational complexity, and this increases 
exponentially with the length of channel response and the number of possible 
states/levels of the modulation employed. The simplified turbo equalisation 
schemes described in Tuchler et al., WO 02/15459 and US 2002/01 10188 consider 

10 systems with a plurality of receive antennas but are incapable of equalisation in 

systems with more than one transmit antenna. It will be appreciated that with 
multiple receive antennas and a single transmit antenna each receive antenna is 
associated with a single channel whereas the situation is much more complex with 
multiple transmit antennas since a combination of signals from different transmit 

1 5 antennas is received. 

There is therefore a need for a reduced computational complexity turbo 
equalization scheme suitable for mitigating both MSI and ISI in a MIMO 
communications system. 

Disclosure of Invention 

20 According to a first aspect of the invention there is therefore provided a 

soft-in-soft-out (SISO) equaliser for use in a receiver of a communications system 
employing a plurality of transmit antennas, the equaliser comprising: at least one 
received signal input for inputting a received signal; a plurality of likelihood value 
inputs, one for each transmit antenna, for inputting a plurality of decoded signal 

25 likelihood values from a SISO decoder; a processor configure to determine from 
the plurality of signal likelihood values an estimated mean and covariance value 
for a signal from each of the transmit antennas; an expected signal determiner 
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coupled to the processor to determine an expected received signal value using the 
mean values; a subtractor coupled to the received signal input to subtract the 
expected received signal value from the received signal to provide a compensated 
signal; a filter coupled to the subtractor to filter the compensated signal to provide 

5 a plurality of estimated transmitted signal values, one for each the transmit 

antenna; a filter coefficient determiner coupled to the processor to determine 
coefficients of the filter using the covariance values; and an output stage coupled 
to the filter to output a plurality of transmitted signal likelihood values, one for 
each the transmit antenna, for input to the SISO decoder. 

10 The equaliser may be used in a receiver of a MIMO communications 

system in which case a plurality of received signal inputs is provided and a MMO 
filter is employed. Broadly speaking use of a multiple output or MIMO linear 
(affine) or transversal filter provides a computational complexity which in 
embodiments increases only linearly with increasing channel impulse response 

1 5 (CIR) length and size of signal constellation. With the aid of soft information 

from the decoder both ISI from the dispersive channels and MSI from the signals 
from the multiple transmit antennas may be jointly mitigated (in embodiments of 
the invention some residual ISI may be present at the equaliser outputs but MSI 
may be substantially removed). Embodiments of the invention thus are able to 

20 . provide a low computational complexity iterative turbo equalisation solution for 
MIMO systems where conventional techniques are impractical Embodiments of 
the equaliser are particularly useful for relatively highly dispersive channels. Any 
suitable conventional decoder may be employed with the equaliser including, but 
not limited to, convolutional, block, turbo, low density parity check or space-time 

25 trellis decoders, or a serially or parallel concatenation of these, depending upon the 
encoder used at the transmitter. 

In embodiments the expected received signal value comprises an ISI 
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(and/or MSI) component of the received signal, for example derived by applying a 
channel response (H) to an expected or estimated transmitted signal (£(X)). 
Preferably the expected received signal value is adjusted by an amount dependent 
upon expected values of estimated transmitted signal values derived from the 

5 signal likelihood values from the SISO decoder, for example expected symbol 
interference values. For an /7-th transmitted symbol interference to the signal Xn 1 
from the i-th transmit antenna may be determined as described later from t*E(x n l ) 
where e 1 is the i-th column of a symbol interference matrix S. 

Preferably the filter coefficients are determined according to a mean square 

10 error cost function, in particular a MMSE function and a filter gain controller is 
provided to control the gain of the filter dependent upon covariance values of 
estimated transmitted signal values, as derived from the signal likelihood values. 
Again the mathematical basis of this is described later. 

To further reduce the computational complexity the filter coefficients may 

1 5 be held constant over a packet or block of received symbols rather than 

recomputed each symbol. Similarly the filter coefficients may be computed in the 
frequency rather than in the time domain, for example in a single-carrier 
modulation system by transforming a packet or block of received symbols to the 
frequency domain for processing in parallel and then afterwards transforming back 

20 into the time domain. Alternatively, in an OFDM (Orthogonal Frequency Division 
Multiplexed) system, where the symbols are in the frequency domain, a 
transformation such as an IFFT (inverse Fast Fourier Transform) is not necessary 
after frequency-domain equalisation. 

A receiver may be provided incorporating both the equaliser and the 

25 decoder and preferably, to improve the independence of the transmitted data, also 
a deinterleaver following the equaliser and an interleaver to interleave the 
likelihood value inputs to the equaliser from the decoder. Such a receiver may be 
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incorporated in, for example, a mobile terminal, base station or access point of a 
communications system with more than one transmit antenna, for example for 
increased data rates and/or, where transmit diversity is employed, for increased 
robustness. 

5 According to another aspect of the invention there is provided a method of 

equalising received data in a communications system employing a receiver with at 
least one receive antenna and a transmitter with a plurality of transmit antennas, 
the method comprising: inputting at least one received signal from the at least one 
receive antenna; inputting a plurality of decoded signal likelihood values, one for 

10 each transmit antenna, from a SISO decoder; determining from the plurality of 

signal likelihood values an estimated mean and covariance value for a signal from 
each of the transmit antennas; determining an expected received signal value using 
the mean values; subtracting the expected received signal value from the received 
signal to provide a compensated signal; determining coefficients for a filter using 

1 5 the covariance values; filtering the compensated signal in accordance with the 

coefficients to provide a plurality of estimated transmitted signal values, one for 
each the transmit antenna; and outputting a plurality of transmitted signal 
likelihood values, one for each the transmit antenna, derived from the estimated 
transmitted signal values, for input to the SISO decoder. 

20 According to a further aspect of the invention there is provided a SISO 

MIMO turbo-equaliser configure for use with a SISO decoder to equalise data 
from a multiple antenna transmitter, the equaliser comprising a multi-dimensional 
transversal filter having a plurality of soft inputs and providing a plurality of soft 
outputs, the equaliser being configure to receive soft information from the SISO 

25 decoder and to use the soft information to adjust coefficients of the transversal 
filter to mitigate MSI and ISL 

In another aspect the invention provides a method of turbo-equalising data 
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10 

in a receiver receiving data from a transmitter configure to transmit data from a 
plurality of transmit antennas simultaneously, the receiver including a SISO 
decoder, the method comprising: inputting a received signal vector comprising a 
block of received signal data including data received from the plurality of transmit 
5 antennas; receiving soft decoded information from the decoder; processing the 

received data using a linear filter having a plurality of soft inputs derived from the 
received signal vector to provide a plurality of soft equaliser outputs, one for each 
transmit antenna, to the SISO decoder; and adjusting coefficients of the linear filter 
in response to the soft decoded information to mitigate MSI and IS! 
10 The invention also provides a method of equalising data in a receiver of a 

communications system with a plurality ni of transmit antennas the method 
comprising: inputting a received signal vector Z„ comprising a block of received 
signal data at an index n\ inputting from a SISO decoder a set of bit likelihood 
values ^(c^ jl / = 1 to ni, j = 1 to m where c^ y , denotes a portion of an m-bit 
1 5 codeword symbol at an index n\ determining expectation and covariance 

cov(x^ , *'„ ) values for estimated transmitted signal values x£ using the values 
^{ c nj)> where x* n denotes an estimated value of a signal transmitted from a 
transmit antenna / at an index w; determining an updated estimated transmitted 
signal value x[ at index n for each transmit antenna i using 
20 xl = K% H [z„ " e ^k))l where E ( z n) = and H is a channel 

impulse response matrix and E(Xr) is an expectation value of an estimated 
transmitted signal vector X„ at index n derived from values E{xl) 9 where e' is the 
ith column of a matrix S and S = h[o„ ;x( „ 2+l _ 1>I; I„ /KWj Q„ iXN]nj f where matrix 
I M is an / x j identity matrix and 0 ixJ is a matrix in which each element is 
25 substantially zero, is the ith column of a filter matrix F„ and F„ is derived from 

cov(x n ,xJ, and where K[ is a scalar filter gain; and outputting equalised 
likelihood values L t {c[ A derived from values x l n . 
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The transmitted symbol index n may index time and/or frequency, 
depending upon the modulation scheme employed. 

The skilled person will recognise that the above-described equalisers and 
methods may be implemented using and/or embodied in processor control code. 

5 Thus in a further aspect the invention provides such code, for example on a carrier 
medium such as a disk, CD- or DVD-ROM, programmed memory such as read- 
only memory (Firmware) or on a data carrier such as an optical or electrical signal 
carrier. Embodiments of the invention may be implemented on a DSP (Digital 
Signal Processor), ASIC (Application Specific Integrated Circuit) or FPGA (Field 

1 0 Programmable Gate Array). Thus the code may comprise conventional program 
code, or micro-code, or, for example, code for setting up or controlling an ASIC or 
FPGA. In some embodiments the code may comprise code for a hardware 
description language such as Verilog (Trade Mark) or VHDL (Very high speed 
integrated circuit Hardware Description Language). As the skilled person will 

1 5 appreciate, processor control code for embodiments of the invention may be 

distributed between a plurality of coupled components in communication with one 
another. 

These and other aspects of the invention will now be further described, by 
way of example only, with reference to the accompanying figures in which: 
20 Brief Description of Drawings 

FIG. 1 shows a known MIMO space-time coded communications system; 
FIG. 2 shows a block diagram of a MIMO channel model; 
FIG. 3 shows a flow diagram of a MMSE MIMO turbo equalisation 
procedure according to an embodiment of the present invention; 
25 FIG. 4 shows a block diagram of a soft-in-soft-out filter-based MEMO 

equalizer according to an embodiment of the present invention; 

FIG. 5 shows a MIMO communications system including a receiver 
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employing the MIMO equaliser of FIG. 4; 

FIG. 6 shows a block diagram of a frequency domain soft-in-soft-out filter- 
based MIMO equalizer according to a further embodiment of the present 
invention; 

5 FIG. 7 shows a receiver incorporating an equaliser configure to operate in 

accordance with an embodiment of an aspect of the present invention; and 

FIG. 8 shows bit error rate performance against received signal to noise 
ratio for a turbo-equalizer and according to an embodiment of the present 
invention, for varying numbers of equalization-decode iterations. 

1 0 Best Mode for Carrying Out the Invention 

Referring to FIG. 2, this shows a block diagram 200 of a MIMO channel 
model. A multi-stream transmitter has first plurality n 1 of transmit antennas 202 
and transmits respective symbols Xn 1 , x n 2 at time n which comprise 

"inputs" to a matrix channel 206. A plurality n 0 of receive antennas 204 provides 

1 5 "outputs" from the « 7 x n 0 MIMO matrix channel in the form of received signals 
Zn ! , Zn 2 . There is a channel with an associated channel response hn^ 

between the ith transmit antenna and jth receive antenna. The received signal at 
each receive antenna also includes a noise component w n J . 

In an n l x n Q MIMO system the MIMO equalizer aims to provide an 

20 estimate of n s transmitted data symbols at every signalling instant. Data from the 
plurality of transmit antennas are transmitted at the same or overlapping times and 
using at the same or overlapping frequencies and thus MSI is introduced, as well 
as ISI from the dispersive wideband channel. In some systems the same or related 
data is transmitted from different transmit antennas in order to provide redundancy 

25 or diversity. In other systems different data streams are transmitted from each 
transmit antenna, for example to provide higher overall data rates. 

Broadly speaking, in embodiments of the invention a turbo equalizer is in 
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r- 

the form of a multi-dimensional transversal filter. Coefficients of the filter are 
adjusted according to a Minimum Mean Square Error (MMSE) criterion to 
mitigate the effects of both Inter-Symbol Interference (ISI) and Multi-Stream 
Interference (MSI) and provide a soft output, with the aid of soft information 

5 relating to multiple transmitted data signals received from a decoder. The 

equalizer processes signals from all the transmit antennas in parallel (although 
potentially in a time-multiplexed manner in processor-based embodiments) is thus 
able to "detect 1 ' and equalise a multistream signal. The soft or likelihood 
information from the decoder is used to determine the covariance and mean of the 

10 multistream transmitted signal for evaluating the coefficients of the equalizer and 
also to determine a mean or expected value of the received signal for MMSE based 
detection. 

In embodiments the SISO equalization may be performed in either or both 
of the time and frequency domains. Initially a time-domain implementation will 
1 5 be described; details of an alternative frequency-domain implementation will be 
given later. 

The elements of the equaliser will generally be implemented by a digital 
signal processor and the structure and operation of a first, time-domain 
embodiment of the equaliser will therefore first be described in mathematical 
20 terms. 

Consider first a signal transmission model for a linear, dispersive and noisy 
communication system with w 7 inputs and n Q outputs (that is a n l xn Q MIMO 
system). The complex-valued baseband equivalent signal model can be 
represented as by equation (1) below: 

25 

L-\ 

Jt=0 

where L is the channel length in symbol periods, where z„ = [z\ z\ • • • z n ° J 
is the signal received by the n Q receive antennas, where 
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10 



15 



20 



25 



w n - \w l n w 2 n — J is the noise component, where 
x n = [ x n x l " • • f * s signal transmitted by the transmit antennas, and 
where n labels a symbol (ie. in this case corresponds to time in symbol periods). 
The dispersive channel matrix is denoted by: 



h" k "* 
h?- 2 



Here, represents the Ath channel tap for a channel link between an /th 
transmit antenna and a jth receive antenna as shown in FIG. 2. 

Consider an embodiment of a transversal filter based equaliser with 
(integer) W coefficients or taps where N = Nl + N2 + 1 (ATI and N2 in effect 
specifying the length of the causal and non-causal portions of the filter). The filter 
employs an observation of a block of received signal of length N in symbol 
periods, which is denoted as: 

Z^HX.+W. (2) 

where 



is the Nn 0 x 1 received signal vector, 
is the Nn Q x 1 received noise vector, 

x„4 



r 



(3) 



(4) 



(5) 



is the (N + L- l)/i y x 1 transmit signal vector which convolves with the 
Nn 0 x(N + L- l)«j channel impulse response (CIR) matrix H given by: 
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15 



20 



25 



H= 



h L-l h L-2 



ho 



(6) 



Determination of the filter coefficients and of a soft transmitted symbol 
estimate by the equaliser will next be considered. 

An optimum linear estimate of the transmitted symbols 
i n = [x\ xl ■■■ x n n ' ^ upon receiving Z„ , based on the second-order statistics 
(that is, means and covariances) has the form given in equation (7) below (see, for 
example, "An Introduction to Signal Detection and Estimation", H. Vincent Poor, 
Second Edition, Springer, ISBN0-387-94173-8 at pages 221 to 233): 

x tt =E{x K ) + A n H [Z a -E{Z n )} (7) 
where, as before, £(•) denotes a mean or expectation value and A„ is the set ofN 
linear filter/equalizer coefficients at symbol (time) n. It can be seen from equation 
(7) that the estimate of the transmitted symbol x n is the mean offset of the filtered 
difference between the received signal vector Z n and the ISI and MSI component 
£(Z n )attimew. 

Minimizing the mean square error cost function £(|x n - i n | 2 ), an optimum 
set of filter coefficients A n is given by: 

A fl = cov(Z w ,Z n )" l cov(Z M ,x w ) (8) 
where the covariance operator cov(x,y) = £(xy H ) - £(x)£(y H ). Equation (7) can be 
expressed in the form: 

x„ ^^(xJ-i-covCx^xJF/K -£(ZJ] (9) 
where the expectation value of the received signal observed at the equalizer taps, 
in terms of the transmitted signal or an estimate of this signal, is given by: 

£(Z>H£(X n ) (9a) 

As will be seen later, £(Z J may be determined from the soft decoder 
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output. In equation (9) the Nn 0 x n. filter matrix F B is defined as 

F = R 77 S 

= |R BW +HR JQr H , 'f S 

where the covariances are defined as 

5 

R zz =cov(Z n) Z„), 

R Jcr =cov(X n) Xj > 

10 IW=cov(W n ,Wj. 

The Nn 0 x w. CIR matrix S introducing the symbol interference from the 
n -th transmitted symbols is given by: 

15 In equation (1 1) matrix I, x , is an i x i identity matrix and matrix 0 inJ 

contains all zeros. 

It can be understood from the form of the equalizer (filter) coefficient 
matrix F„ expressed in equation (10), which includes the channel matrix Has 
defined in equation (6), that the equalizer not only performs equalization on the ISI 

20 but also mitigates MSI. 

The mean and covariance of the transmitted signal, E(x n ) andcov(x /J ,x rt ), 
which are used to compute the filter coefficients and the estimated transmitted 
signal, may be obtained using equations (16) and (17) below. The block-diagonal 
covariance matrix comprises (N+L-I) block-diagonal blocks of cov(x„,x n ) 

25 as follows: 

C0V ( X n-N7-U\ > X n-N2-L+\ ) 

The CIR matrix H may be obtained from a MIMO channel estimation 
block in the receiver in a conventional manner. Thus in one embodiment known 
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pilot or training sequences are periodically inserted into the signal from each 
transmit antenna and at the receiver these known sequences are encoded and 
provided to a channel estimator together with one (or more) input signal streams 
from the one (or more) receive antennas. Providing the training sequences in the 

5 signal streams have a low cross-correlation coefficient the channel estimator is 
then able to determine the CIR matrix H by, in effect, correlating the encoded 
known data sequence for each stream with the actual received data. An example 
of such a channel estimator is described in Ye Geoffrey Li, "Simplified channel 
estimation for OFDM systems with multiple transmit antennas", IEEE 

10 Transactions on Wireless Communications, Vol. 1, No. 1, pg. 67, Jan 2002, which 
is hereby incorporated by reference. 

Advantageously the channel estimation may also be performed iteratively 
by using the estimated transmitted symbols x' n as additional training symbols. For 
background information relating to this reference may be made to Tetsushi Abe 

1 5 and Tad Matsumoto, "Space-Time Turbo Equalization and Symbol Detection in 
Frequency Selective MIMO Channels" in: Proc. Veh. Techn. Conference, IEEE 
VTS 5 th Vol. 2. pg 1230-1234, 2001, also hereby incorporated by reference. Thus, 
for example, multiple, deinterleaved outputs from the decoder (one for each 
transmit antenna) may be re-encoded and provided to the channel estimator 

20 similarly to a known training sequence. In this way, an updated estimated value of 
H may be obtained on every turbo equalization iteration and this may in turn be 
used to compute the equalizer coefficients. 

For good convergence behaviour the soft transmitted symbol information 
K provided by the equalizer should be independent from the soft transmitted 

25 symbol information from the decoder in the iterative structure. At the transmitter a 

sequence of m bits, where m is the number of bits per symbol, is encoded as a 
codeword e {+1,-1}, j = 1,— ,jw for transmit antenna The decoder provides 
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soft log likelihood ratio information L{c' nj \j = !,-••, m where 

L { C nj) ~ lli { P ( C nJ = + l ) /P { C n.j = ~ ! )} Wflich > fr0IT1 ^ P° illt ° f vieW ° f the 

equaliser, may be considered as a priori information. To provide good 
convergence behaviour the likelihood information for bits c^ J9 j = 1,.--,to 

5 (corresponding to symbol x[ ) which are available to the equalizer is set to 

Z,(cj, j y qualizer = o, y = 1, • • • , m to derive an equation for x* n . (During the turbo 
equalisation these values may also be initialised to zero). Then, eliminating the 
dependence of the equalisation on L^l ^j = 1,—,/w provided by the decoder, and 
setting the expectation E(x i n y quabzer =o m ^ t h e covariance cov^,*^)^*'*" = 5] 4 , 

10 where S 2 x i is an initial value, for example unity (depending upon the relative gain 
of the transmit antennas), for the signal variance of transmitted symbols from 
transmit antenna i , the i th component of the estimated transmitted symbol i„ in 
equation (9) can be expressed as scaled version of F„ (using the matrix inversion 
lemma), as shown in equation (12) below: 

15 

^=<fr[z n -(£(Z M )-e'£(x:))] 02) 
where the filter gain K l n is given by: 

^^^fe-covfe.x^e')" (13) 
20 In equations (12) and (13) is the / th column of the filter matrix 

F n defined in equation (10) and e' is the /th column of the symbol interference 
matrix S given in equation (11). 

Next a reduced complexity method of determining R^" 1 is described, for 
computing the filter matrix F n via equation (10). In order to avoid performing a 
25 direct matrix inversion of the Nn 0 x Nn Q covariance matrix R ffi to compute the 

filter matrix F B for every estimate of symbols i n , a time-recursive update 
algorithm may be employed. The algorithm is similar to one described in the 
Tuchler et al. paper (ibid) and exploits the existence of common submatrices 
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within a partitioned covariance matrix. The matrix R^ 1 at a time n+l may be 
determined by the manipulation of submatrices as set out in equation (14): 



(14) 



These submatrices may in turn be determined based upon a previous 
at a time n , as shown by equations (15) below: 



-i 



- H 
Up Up 



05) 



The sizes of the matrices © w , cr„ , Up , , are {N - \)n Q x n Q , 
10 n 0 xn 0> n 0 xn 0 , (N- \)n Q x w 0 and (tf - l)w 0 x (TV - l)w 0 > respectively. It can 

therefore be seen that only the matrix inversion of n Q xn Q matrices u p and 
( <j N - g> N H tu N ) is required to update the matrix R^" 1 . The term R^ is the noise 
covariance and a value for this may be determined, for example, theoretically from 
the receiver bandwidth (particularly where the front end filter has a relatively 
1 5 sharp cut-off; see also 3 GPP TS25.2 1 5 v5.2. 1 for background on received power), 

or by a measurement of the level of noise (and/or interference) at the receiver, or 
by a combination of both these techniques. Where the noise levels at the receive 
antennas are similar R^ may approximate to a fraction of I. 

The soft (likelihood) value inputs to the equaliser from the decoder will 
20 now be discussed. 
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The soft information provided by the decoder in terms of the likelihood 
values are utilized to provide the mean E{x' n ) and covariance cov^xjjof the 
transmitted symbols, which are required to compute x' n in Equation (9). The mean 
and covariance values are obtained as shown below: 

£fo)=Z<v(<=«*) 06) 



a k *S 



and 



10 At the transmitter a sequence of m bits c is mapped to a symbol x selected 

from a Q - 2 m -ary symbol alphabet S = \p u a 2i '- 7 a Q } in which each symbol a* 
corresponds to an m bit pattern = [s k ] s k 2 J^J, e {+l,-l}. The 
probabilities of the transmitted symbols are given by: 

4l=aJ=n; =1 p(< /= 5j 

15 =n; =l l/2 ( f + 5 U tanh(i(< y )/2)} (18) 

1=1, — = 1,— ,6 

The soft (likelihood) value outputs from the equaliser to the decoder will 
next be considered. 

At the output of the equalizer, it is assumed the probability density 
functions of the estimated symbols p{x i n | c' n = s*) are Gaussian with mean j^ nJ( 
20 and variance <j 2 J . The log-likelihood ratios of the bits contained in the transmitted 
symbols are then given by: 
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10 



15 



4(<>ln 



£ ex p 



-In 



I ex P~ rr— a - + Aj 



'n,k 



)) 



(19) 



v/VW 

Based on the derivation obtained for the estimate of the transmitted symbol 
x' n in equation (12) the mean and variance of the equalizer output may then be 
written: 



M« =*i/;" Wn \< =a k )-(E{Z n )-* i E(x[ 



(20) 



and 



(21) 



cov(Z,„ZJ<=a t )fi 

= K? (ffRJl - cov $ t'SteX) 

The logarithmic summation of exponentials used in equation (19) to 
compute the equalizer's soft output, ln^exp(A t )j, may advantageously be 
computed using the Jacobian logarithmic relationship: 



f n 



\k=\ J 

yU,^)«maxU J A 2 )+/ c ^-^|) 

Here / c (x) is a logarithmic correction function that can be tabulated in a 
look-up table such as that given in Table 1 below. 
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Table 1 



X 


/.W 


x>3.7 


0.0 


3.7 £x> 2.25 


0.05 


2.25 >x> 1.5 


0.15 


1.5>x>1.05 


0.25 


1.05 £ x> 0.7 


0.35 


0.7 £x> 0.43 


0.45 


0.43 >x> 0.2 


0.55 


0.2 >x 


0.65 



In this way the computational complexity of technique may be reduced by 
replacing the use of logarithmic and exponential functions in the likelihood 
calculation with a simple look-up table and "maxO" function. 
5 To further reduce the computation complexity, the equalizer coefficients, 

that is terms F n and K[ defined in equations (10) and (13) respectively, may be 
set to be non- varying with time n. An average of the covariance matrix 

=-^Z C0V ( X "> X «) and cov (<> x J = ^Z cov fc><) overareceived 
symbol packet of P symbols, may instead be used to compute the coefficients in 
10 . Equation (10) and (13). These time invariant equalizer coefficients F and K* may 
then be used to facilitate a reduced complexity computation of the transmitted 
symbol x[ in Equation (12). 

FIG. 3 shows a flow diagram of an implementation the above-described 
turbo equalisation procedure. System parameters are determined or chosen at step 
15 S800 including, for example, the symbol alphabet S and the number of filter 

coefficients N=Nl+N2(N\ and N2 are typically chosen such that N\=N2=L). 
Strictly speaking, however, this is not part of the actual equalisation procedure. 

At step S802 equaliser parameters are initialised, for example by setting 
E{x i n ) =0, cov^x) and L(<J=0 (a log likelihood ratio 1^ of 0 

20 implying c=+l and c=-l are equiprobable), then at step S804 a block of received 
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signal Z n is input from one or more receivers and, at step S806, a estimated OR 
matrix H is input from a channel estimator. 

At step S808 the mean E(x i n ) and covariance cov^.xjjof the transmitted 
symbols are obtained, according (after the first iteration) to the a priori 

5 information of the transmitted bits from the SISO decoder 

L(cl = 1,-- -jTijJ = !,*••,/» (where m is the number of bits per symbol) using 
equations (16), (17) and (18). Then, at step S812, using the mean E{x i n ) and 
covariance cov^^xj,) values of the transmitted symbols the equalizer coefficients 
F n and K[ are obtained (using equations (10) and (13) respectively) and the mean 

10 of the received signal E(Z n ) is determined using equation (9a). The estimated 
transmitted signal i„ is then given by equation (12), and the "extrinsic" 
information on the transmitted symbols, L g [c' nJ ) , which will provide a soft input 
to the SISO decoder from the equaliser, is obtained using equation (19) given the 
estimated transmitted symbols and the a priori information relating to the 

1 5 transmitted bits. 

Deinterleaving is performed at step S814, corresponding to the interleaving 
performed after encoding and before transmission of the data at the transmitter. 
Thus the extrinsic information from the equalizer is deinterleaved before passing 
to the decoder. Then at step S816 the extrinsic information relating to the 

20 transmitted bits is used as the a priori knowledge at the decoding stage to provide 
the decoder's extrinsic information. The deinterleaved signals are decoded in 
accordance with the encoding used at the transmitter. A check is then made, at 
step S818, to determine whether a termination criterion has been reached and, if it 
has not, the procedure loops back to reinterleave the extrinsic information from the 

25 decoder and to pass this back to the equalizer via step S808. The termination 

criterion may comprise, for example, a predetermined number of iterations or a 
determination of whether the decoder soft output, that indicates the reliability of 
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the decoded data, is more than a threshold value. If the termination criterion has 
been met the equalized and decoded bits at the final iteration are provided as an 
output estimate of the transmitted data bits. 

FIG. 4 shows a block diagram of a soft-in-soft-out filter-based MMO 

5 equalizer 400 configure to operate according to the procedure of FIG. 3. 

Referring to FIG. 4 the equaliser comprises an input 402 to receive a 
plurality of received signal blocks from a corresponding plurality of receive 
antennas, these received signal inputs providing a first set of inputs to a subtractor 
404. A second set of inputs to the equaliser 400 comprises a set of soft or 

10 likelihood values 406 from a soft-in-soft-out decoder, to provide a priori 

transmitted bit information to the equaliser. This information is processed by 
block 405 to determine e(x[) and cov(x^,xj,) and these (estimated) mean (or 
expectation) and covariance values are used by block 408 to perform a calculation 
[E(Z H )-e i E{x i Jj to determine an expected set of received signal values. This 

15 provides a second set of inputs 410 to subtractor 404 for subtraction from the first 

set of inputs 402 to provide a "compensated" signal 412 to a linear (affine) or 
transversal filter 414. The soft likelihood values 406 from SISO decoder are used 
by a filter coefficient calculation block 416 to determine a set of coefficients for 
filter 414. A MIMO channel estimator 418 provides a channel estimate input 420 

20 to this filter coefficient calculation block 416. The output of filter 414 comprises a 
set of estimated transmitted signal values 422, one for each transmit antenna, 
which are in turn provided to a soft decision block 424 which provides a 
corresponding plurality of outputs 426 comprising soft (or likelihood) transmitted 
bit values which are suitable for providing a soft input to the SISO decoder. 

25 FIG. 5 shows a MIMO communication system 500 including a MEMO 

transmitter 502 and a MIMO receiver 504 communicating via a MIMO channel 
506. The receiver 504 incorporates a MIMO SISO (MMSE) equaliser 400, for 
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example as described with reference to FIG. 4. The transmitter 502 has a data 
input 508 providing an input to a space-time and/or channel encoder 510 (in other 
embodiments other coding schemes may be employed). The encoder 510 provides 
a plurality of outputs 5 12 to an interleaver 514 which, in turn, provides signals to a 

5 plurality of rf output stages and thence to a corresponding plurality of transmit 

antennas. In FIG. 5, for clarity these rf stages and transmit antennas are not shown 
explicitly but instead a plurality of outputs 516 comprising the plurality of 
transmitted signals is shown providing an output from interleaver 514 which, via 
MIMO channel 506, provides a (different) plurality of inputs to receiver 504. It 

10 will be recognised that the communication system of FIG. 5 employs a single 

transmitter to provide a plurality of transmit output streams, either for redundancy 
or increased bit rate. The transmissions from the plurality of transmit antennas 
may, for example, share the frequency or overlap in frequency and/or overlap in 
time. This is different to a communication system employing a plurality of users 

1 5 with frequency and/or time domain controlled access in which, generally speaking, 
it is preferred to assign different frequencies and/or time slots to different uses. 

Continuing to refer to FIG. 5, a plurality of receive antennas coupled to a 
corresponding plurality of rf receiver front ends (not shown for clarity in FIG. 5) 
provides a plurality of inputs to MIMO SISO MMSE equaliser 400. The soft 

20 output 426 from the equaliser 400 are deinterleaved by a deinterleaver 5 1 8. and 
then provided to a space-time/channel decoder 520. The decoder accepts a 
plurality of inputs, one for each signal stream from a transmit antenna, and 
provides a corresponding plurality of outputs 522a, b which are either provided to 
a bit interleaver 524 and returned to equaliser 400 for a further equalisation- 

25 decoding iteration or, if a termination criterion has been reached, output as 

estimated data. It will be appreciated that equaliser needs complete transmitted 
symbols from decoder 520, that is where, for example, error check bit such as 
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parity bits have been included transmit symbols including these parity bits should 
be provided to the equaliser. The MEMO channel estimator 418 of FIG. 4 (not 
shown in FIG. 5) may receive a set of inputs from bit interleaver 524 and output an 
estimate of H to equaliser 400. 

5 Referring now to FIG 6 this shows a block diagram of a soft-in-soft-out 

filter-based MIMO equalizer 600 configure to operate in the frequency domain. 
The main elements of equaliser 600 are similar to those of FIG. 4 and are indicated 
by like reference numerals. The main differences are that calculation blocks 408 
and 416 perform frequency domain calculations and that subtracter 404 and filter 

10 414 operate in the frequency rather than the time domain (though, for simplicity of 

comparison, the same reference numerals have been employed and the channel 
estimator has been omitted). Thus the main additions are a serial-to-parallel 
conversion block 602 and fast Fourier transform (FFT) block 604 prior to 
subtracter 404 and filter 414, and an inverse fast Fourier transform (FFT 1 ) 606 

15 and a parallel-to-serial conversion block 608 following filter 414. These blocks 

convert the received signals from the time to frequency domain for processing and 
back to the time domain after processing. 

Background information relating to SISO equalisation in the frequency 
domain can be found in Tuchler, M., Hagenauer, J., "Turbo equalization' using 

20 frequency domain equalizers" in: Proa Allerton Conference, Monticello, EL, USA, 

Oct 2000, to which reference may be made. It will be appreciated, however, that 
this reference only discloses frequency domain turbo equalisation for a single- 
transmit-antenna-single-receive-antenna system. Other background information 
can be found in M V. Clark, "Adaptive Frequency-Domain Equalization and 

25 Diversity Combining for Broadband Wireless Communications", IEEE JSAC, vol. 

16, no. 8, pp. 1385-1395, Oct. 1998, and D. Falconer, S. L. Ariyavistakul, A. 
Benyamin-Seeyar, B. Eidson, "Frequency Domain Equalization for Single-Carrier 
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Broadband Wireless Systems", IEEE Communications Magazine, pp. 58-66, 
April 2002. 

Continuing to refer to the example of a frequency domain implementation 
of the equaliser shown in FIG. 6, the signals from the receive antenna(s) are 

5 transformed into the frequency domain by Fast Fourier Transform (FFT) on a per 

FTT block/packet basis and the spectra of the received signals are equalized by 
multiplying the spectrum of each branch with the frequency domain coefficients of 
the equalizer. Frequency domain equalization again reduces the computational 
complexity of the equalization task by processing an entire FFT block/packet of 

10 received symbols at once in the frequency domain instead of symbol-by-symbol in 
the time domain although it will be appreciated that because of this the equalizer 
coefficients are time-invariant per FFT block. Furthermore a cyclic extension (for 
example, prefix) should be added to the transmitted symbol block/packet before 
transmission (to provide an integral number of cycles within the FFT interval) 

15 although this is discarded at the receiver. The frequency domain implementation 

of the turbo equaliser is suitable for both time and/or frequency domain coded data, 
and may be used, for example, with MIMO OFDM transmissions. 

Referring now to FIG. 7, this shows an embodiment of a receiver 700 
incorporating a equaliser configured to operate as described above, in particular to 

20 operate as described with reference to FIGS. 3 and 4. 

Receiver 700 comprises one or more receive antennas 702a, b (of which 
two are shown in the illustrated embodiment) each coupled to a respective rf front 
end 704a, b, for example similar to the rf front end of FIG. 2a, and thence to a 
respective analogue-to-digital converter 706a,b and to a digital signal processor 

25 (DSP) 708. DSP 708 will typically include one or more processors 708a (for 

example, for a parallel implementation of filter 414) and some working memory 
708b. The DSP 708 has a data output 710 and an address, data and control bus 
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712 to couple the DSP to permanent program memory 714 such as flash RAM or 
ROM. Permanent program memory 714 stores code and optionally data structures 
or data structure definitions for DSP 708. 

As illustrated program memory 714 includes SISO equalisation code 714a 

5 comprising [E(Z n )-e i E(x i Jj calculation code, subtraction code, filter coefficient 
calculation code, linear (transversal) filter code and soft decision output code to, 
when running on DSP 708, implement the corresponding functions as described in 
detail above. Program memory 714 also includes MIMO channel estimation code 
714b to provide a MIMO CIR estimate H, de-interleaver code 714c, interleaver 

10 code 714d, and space-time (or other) SISO decoder code 714e. Examples of 

suitable code 714e are well known to those skilled in the art. Optionally the code 
in permanent program memory 714 may be provided on a carrier such as an 
optical or electrical signal carrier or, as illustrated in FIG. 7, a floppy disk 716. 
The data output 710 from DSP 708 is provided to further data processing elements 

15 of receiver 700 (not shown in FIG. 7) as desired. These may include a block error 
decoder such as a Reed- Solomon decoder (although this could be part of the turbo 
decoder), and a baseband data processor for implementing higher level protocols. 

The receiver front-end will generally be implemented in hardware whilst 
the receiver processing will usually be implemented at least partially in software 

20 although one or more ASICs and/or FPGAs may also be employed. The skilled 
person will recognise that all the functions of the receiver could be performed in 
hardware and that the exact point at which the signal is digitised in a software 
radio will generally depend upon a cost/complexity/power consumption trade-off. 
FIG. 8 shows bit error rate (BER) performance against received signal to 

25 noise ratio (Eb/No) for an embodiment of the turbo-equalizer, for varying numbers 
of equalization-decode iterations. In FIG. 8 curve 800 shows a theoretical result 
for a single transmit single receive antenna system with a non dispersive channel 
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and the remaining curves show the 8PSK performance of an MMSE-based MIMO 
SISO turbo-equalizer for a two-transmit-antenna transmitter, two-receive-antenna 
receiver MEMO system with a dispersive channel of 

h ] 0 A = 0.707, ti? = 0.707, /£ 2 = 0.707, =0.707, 
*f 1 = 0.707,/? 2 ' 1 = 0.707.A* 2 = 0.9579,// 2 * 2 = 0.2874 " 

The system utilizes half-rate (5,7) convolutional coding and the equalizer 

parameters are Nl=9 and Af2=5. More specifically curves 802, 804, 806, 808 and 

810 show, respectively, curves for one, two, three, four and five turbo equalisation 

iterations. It can be seen that performance increases significantly with the number 

10 of iterations and approaches that of a non-dispersive channel, although with large 

numbers of iterations the incremental performance increase is reduced. The 
reduction in incremental increase can be used to determine a termination criterion 
for a potential application, for example limiting to 2, 5 or 10 iterations depending, 
for example, upon speed-processing power trade-offs. 

1 5 Applications of the invention have been mainly described in the context of 

a MIMO system with time domain coding but embodiments of the invention are 
also useful in frequency domain coded systems such as MEMO-OFDM . 
(Orthogonal Frequency Division Multiplexed) systems. For example the invention 
may be employed with the European Hiperlan/2 or US IEEE 802.11a standards for 

20 54Mbps wireless networks. Embodiments of the invention may also be employed 
in non-wireless applications such as magnetic or optical disk drive read head 
circuitry where, for example, multiple layers of a disk in effect act as multiple 
transmitters, one or more heads receiving read data influenced by "transmitted" 
signals from more than one layer. 

25 No doubt many other effective alternatives will occur to the skilled person. 

It will be understood that the invention is not limited to the described embodiments 
and encompasses modifications apparent to those skilled in the art lying within the 
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spirit and scope of the claims appended hereto. 
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CLAIMS 

1 . A soft-in-soft-out (SISO) equaliser for use in a receiver of a 
communications system employing a plurality of transmit antennas, the equaliser 
comprising: 

5 at least one received signal input for inputting a received signal; 

a plurality of likelihood value inputs, one for each transmit antenna, for 
inputting a plurality of decoded signal likelihood values from a SISO decoder; 

a processor configure to determine from said plurality of signal likelihood 
values an estimated mean and covariance value for a signal from each of said 
10 transmit antennas; 

an expected signal determiner coupled to said processor to determine an 
expected received signal value using said mean values; 

a subtractor coupled to said received signal input to subtract said expected 
received signal value from said received signal to provide a compensated signal; 
15 a filter coupled to said subtractor to filter said compensated signal to 

1 provide a plurality of estimated transmitted signal values, one for each said 
transmit antenna; 

a filter coefficient determiner coupled to said processor to determine 
coefficients of said filter using said covariance values; and 
20 an output stage coupled to said filter to output a plurality of transmitted 

signal likelihood values, one for each said transmit antenna, for input to said SISO 
decoder. 

2. A SISO equaliser as claimed in claim 1 further comprising an adjuster 
to adjust said expected received signal value by an amount dependent upon said 

25 estimated mean values. 

3. A SISO equaliser as claimed in claim 2 further comprising a filter gain 
controller configure to control the gain of said filter dependent upon said estimated 
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covariance values. 

4. A SISO equaliser as claimed in any preceding claim wherein said filter 
comprises a linear or transversal filter. 

5. A SISO equaliser as claimed in any preceding claim wherein said filter 
5 coefficient determiner is configure to determine said filter coefficients according 

to a mean square error cost function. 

6. A SISO equaliser as claimed in claim 1 wherein said filter coefficient 
determiner is configure to determine said filter coefficients responsive to 
covariance values of estimated transmitted signal values derived from said signal 

1 0 likelihood values from said SISO decoder. 

7. A SISO equaliser as claimed in any preceding claim configure to utilise 
substantially constant filter coefficients for equalising a block or packet of 
received data symbols. 

8. A SISO equaliser as claimed in claim 7 wherein said filter coefficient 
1 5 determiner is configure to operate in the frequency domain, said equaliser further 

comprising Fourier transform means prior to said filter and inverse Fourier 
transform means following said filter. 

9. A SISO equaliser as claimed in any preceding claim for use in a 
receiver for receiving data from a transmitter configure to transmit data from a 

20 plurality of transmit antennas simultaneously, the equaliser comprising a plurality 
of received signal inputs each coupled to said subtractor for inputting a plurality of 
received signals one from each of a plurality of receive antennas, and wherein said 
subtractor is configure to subtract a plurality of expected signal values relating to 
said plurality of receive antennas from said plurality of received signals to provide 

25 compensated signals to said filter. 

10. A receiver incorporating the SISO equaliser and decoder of any 
preceding claim. 
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1 1 . A SISO MEMO turbo-equaliser configure for use with a SISO decoder 
to equalise data from a multiple antenna transmitter, the equaliser comprising a 
multi-dimensional transversal filter having a plurality of soft inputs and providing 
a plurality of soft outputs, the equaliser being configure to receive soft information 

5 from said SISO decoder and to use said soft information to adjust coefficients of 
said transversal filter to mitigate MSI and ISL 

12. A SISO MMO turbo-equaliser as claimed in claim 1 1 further 
configure to adjust said transversal filter coefficients in accordance with a 
minimum mean square (MMSE) criterion. 

10 13. A SISO MIMO turbo-equaliser as claimed in claim 1 1 or 12 wherein 

said filter operates in the frequency domain, and said coefficients comprise 
frequency domain coefficients, the equaliser further comprising Fourier transform 
means preceding said filter and inverse Fourier transform means following said 
filter. 

1 5 14. A SISO MIMO turbo-equaliser as claimed in claim 1 1 or 12 wherein 

said filter operates in the time domain and wherein said coefficients comprise 
coefficients which are substantially time invariant over a symbol packet 
comprising a plurality of received symbols. 

15. A receiver including the SISO MIMO turbo-equaliser of any one of 
20 claims 1 1 to 14 and said SISO decoder, said decoder being configure to receive 

said soft equaliser outputs as inputs. 

16. A method of equalising received data in a communications system 
employing a receiver with at least one receive antenna and a transmitter with a 
plurality of transmit antennas, the method comprising: 

25 inputting at least one received signal from said at least one receive antenna; 

inputting a plurality of decoded signal likelihood values, one for each 
transmit antenna, from a SISO decoder; 
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determining from said plurality of signal likelihood values an estimated 
mean and covariance value for a signal from each of said transmit antennas; 

determining an expected received signal value using said mean values; 

subtracting said expected received signal value from said received signal to 
5 provide a compensated signal; 

determining coefficients for a filter using said covariance values; 

filtering said compensated signal in accordance with said coefficients to 
provide a plurality of estimated transmitted signal values, one for each said 
transmit antenna; and 

10 outputting a plurality of transmitted signal likelihood values, one for each 

said transmit antenna, derived from said estimated transmitted signal values, for 
input to said SISO decoder. 

17. A method as claimed in claim 16 further comprising adjusting said 
expected received signal value by an amount dependent upon said estimated mean 

15 values. 

18. A method as claimed in claim 17 further comprising controlling the 
gain of said filter dependent upon said estimated covariance values. 

19. A method as claimed in claim 16, 17 or 18 wherein said filtering 
comprises linear or transversal filtering. 

20 20. A method as claimed in any one of claims 16 to 19 wherein said 

determining substantially minimises a mean square estimated transmitted signal 
error cost function. 

21. A method as claimed in claim 16 wherein said determining is 
responsive to covariance values of estimated transmitted signal values derived 

25 from said signal likelihood values from said SISO decoder. 

22. A method as claimed in any one of claims 16 to 21 further comprising 
maintaining said filter coefficients substantially constant over a block or packet of 
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received data symbols. 

23. A method as claimed in any one of claims 16 to 22 wherein said 
filtering is performed in the frequency domain. 

24. A method as claimed in claim 23 further comprising transforming said 
5 received signal into the frequency domain prior to said filtering and transforming 

the result of said filtering into the time domain following said filtering. 

25. A method as claimed in any one of claim 16 to 24 for use in a receiver 
receiving data from a transmitter configure to transmit data from a plurality of 
transmit antennas simultaneously, the method comprising inputting a plurality of 

10 received signals from a plurality of receive antennas, and wherein said subtracting 
comprises subtracting a plurality of expected signal values relating to said plurality 
of receive antennas from said plurality of received signals to provide compensated 
signals for said filtering. 

26. A method of turbo equalisation comprising iteratively equalising data 
1 5 according to any one of claims 1 8 to 25 and decoding said transmitted signal 

likelihood values from said equalising to provide decoded signal likelihood values 
for said equalising. 

27. A method of turbo-equalising data in a receiver receiving data from a 
transmitter configure to transmit data from a plurality of transmit antennas 

20 simultaneously, the receiver including a SISO decoder, the method comprising: 
inputting a received signal vector comprising a block of received signal 
data including data received from said plurality of transmit antennas; 
receiving soft decoded information from said decoder; 
processing said received data using a linear filter having a plurality of soft 
25 inputs derived from said received signal vector to provide a plurality of soft 

equaliser outputs, one for each transmit antenna, to said SISO decoder; and 

adjusting coefficients of said linear filter in response to said soft decoded 
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information to mitigate MSI and ISI. 

28. A method of equalising data in a receiver of a communications system 
with a plurality ni of transmit antennas the method comprising: 

inputting a received signal vector Z„ comprising a block of received signal 
5 data at an index n\ 

inputting from a SISO decoder a set of bit likelihood values £(cj t/ )i i - 1 
to nj, j = 1 to m where c* n y , denotes a portion of an m-bit codeword symbol at an 
index n\ 

determining expectation E(xl) and co variance cov(x^x^) values for 
10 estimated transmitted signal values x[ using the values l(cJ, s \ where x? n denotes 

an estimated value of a signal transmitted from a transmit antenna / at an index n\ 

determining an updated estimated transmitted signal value x x n at index n 
for each transmit antenna / using x[ = K 9 n t § * [z„ - ^(zj-e'^x*))] where 
E(Z n ) = H£(X n ) and H is a channel impulse response matrix and E(S.r) is an 
1 5 expectation value of an estimated transmitted signal vector X„ at index n derived 

from values E[x* n ) y where e' is the ith column of a matrix S and 

f where matrix l M is an i xi identity matrix 
and 0 ix/ is a matrix in which each element is substantially zero, f l n is the ith 
column of a filter matrix F„ and F„ is derived from cov(x n ,x„), and where K* H is a 
20 scalar filter gain; and 

outputting equalised likelihood values L t derived from values x* n . 

29. A method of equalising data as claimed in claim 28 further comprising 
decoding said values t0 provide said values iiflj). 

30. A method as claimed in claim 29 comprising iteratively equalising and 
25 decoding to determine values for L 9 (cj,,.) and Z,(c^ ) until a termination criterion 

is reached. 

31. A method as claimed in claim 30 further comprising deinterleaving 
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said values and interleaving said l(c^ ; ) values. 

32. A method as claimed in any one of claims 28 to 3 1 wherein F„ is 
derived from [R W w + HR^H^S where Ra^= cov(x„,x n ) and R W w is a noise 
covariance matrix. 

5 33. A method as claimed in any one of claims 28 to 32 wherein n indexes 

time. 

34. A method as claimed in any one of claims 28 to 32 wherein n indexes 
frequency. 

35. Processor control code to, when running, implement the method of any 
1 0 one of claims 1 6 to 34. 

36. A carrier carrying the processor control code of claim 35. 

37. An equaliser configure to operate in accordance with the method of 
any one of claims 16 to 34. 

38. A data communications receiver including the equaliser of claim 37. 
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